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ABSTRACT 

The nearby A4-type star Fomalhaut hosts a debris belt in the form of an eccentric ring, which is thought to 
be caused by dynamical influence from a giant planet companion. In 2008, a detection of a point-source 
inside the inner edge of the ring was reported and was interpreted as a direct image of the planet, named 
Fomalhaut b. The detection was made at ^600-800 nm, but no corresponding signatures were found in 
the near-infrared range, where the bulk emission of such a planet should be expected. Here we present 
deep observations of Fomalhaut with SpitzerflRAC at 4.5 /im, using a novel PSF subtraction technique 
based on ADI and LOCI, in order to substantially improve the Spitzer contrast at small separations. The 
results provide more than an order of magnitude improvement in the upper flux limit of Fomalhaut b and 
exclude the possibility that any flux from a giant planet surface contributes to the observed flux at visible 
wavelengths. This renders any direct connection between the observed light source and the dynamically 
inferred giant planet highly unlikely. We discuss several possible interpretations of the total body of 
observations of the Fomalhaut system, and find that the interpretation that best matches the available data 
for the observed source is scattered light from transient or semi-transient dust cloud. 

Subject headings: circumstellar matter — planetary systems — stars: early-type 



1. Introduction 

Direct imaging of exoplanets has been a field in 
rapid development over the pas t few years, with de- 



could potentially be a p lanet in the process of forming 
dKraus & Irelandll201ll) . 



tections of several planets (e.g. Lagrange et al. 2010; 



Marois et al. 2010) a nd low-mass brown-d warfs (e.g. 
Chauvinetal.ll2005t iThalmann etai]|2009h that have 
been enabled by the advent of high-quality adaptive 
optics correction and sophisticated PSF subtraction 
techniques. This has allowed for v arious kinds of char- 
acterization of such systems (e.g . Jjanson et alj| 2010; 
Bowler et al. 2010t Currie et al. 201 lb Ijanson et al.l 
201lb . and opened up new domains for the types of 
planets that can be studied; for instance, observations 
within the disk gap of the transitional disk LkCal5 



Among these discoveries, one claimed planet de- 
tection that stands out as peculiar in many ways is 
that of Fomalhaut b dKalas et al.ll2008i hereafter K08). 
The presence of a planet around Fomalhaut has been 
predicted on the bas is of the geometr y of the debris 
disk in the system ( Kalas et al. 2005), which has a 
sharp inner edge and a center that is offset from that 
of the star. This implies that it must be eccentric, 
which could in turn be indicative of the presence of 
a planetary companion (or several) e xerting a grav- 
itation al influence on the disk (e.g. lOuillenl [2006; 
Chiang et al.l 12009b. although alternative interpreta- 



(e.g. Espaillat et alJl2007 : Thalmann et al]l201(3) using 
sparse aperture masking has recently revealed what 
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tions do exist (Jalali & Tremainell201 1 ). Hence, when 
a point-source was discovered in two epochs within the 
disk gap (K08), with a direction of motion largely par- 
allel to the disk edge, it was assumed that this was an 
image of the predicted disk-perturbing planet. How- 
ever, the observed properties of the point-source are 
hard to consolidate with such an interpretation. Un- 
like the other detections mentioned above, which were 
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made at near-infrared wavelengths where young sub- 
stellar objects radiate the bulk of their energy, the Fo- 
malhaut b candidate is detected only at visible wave- 
lengths, where the expected emission is near zero. 
No corresponding near-infrared radiation has so far 



been detected, despite sev eral attempts dKalas et al 
120081 iMarengo etal]l2009h . Several alternative inter- 
pretations of the observed properties have been made, 
which will be discussed in Sect. |4] Regardless of 
interpretation, it remains the case that the best way 
to increase our understanding of the system and test 
whether a planet is associated with the observed point- 
source is to better constrain its near-infrared proper- 
ties. 

Motivated by this, we have performed a study with 
the Spitzer Space Telescope in order to improve the 
detection limits at 4.5 p,m, which is the wavelength 
range where Fomalhaut b is expected to emit its peak 
flux. We describe our observational methods and data 
reduction in Sect. [2] our various data analysis ap- 
proaches and results in Sect. [3] and discuss the im- 
plications for the Fomalhaut system in Sect. |4] 

2. Observations and Data Reduction 

Our observations were taken with the Infrared Ar- 



ray Camera (IRAC; iFazio et al.ll2004 of the Spitzer 
Space Telescope as part of program 70009, and con- 
sist of eight individual runs spread over cycle 7, from 
August 2010 through July 2011 (see Table QJ- Each 
run consists of 48 exposures, structured as a cycle of a 
12-point Reuleaux dither pattern with four exposures 
per position, which enables efficient spatial oversam- 
pling and bad pixel removal. The individual exposures 
have integration times of 10.4 s each with an execution 
time of 12 s, leaving the primary saturated in individ- 
ual frames. All observations were taken in the 4.5 /im 
band where the peak flux of Fomalhaut b is expected. 
The observing strategy was optimized for angular dif- 
ferential imaging (ADI) purposes, with a large spread 
in telescope roll angles. Since active rotation around 
the optical axis of the telescope is not possible, we ex- 
ploited the fact that nominal rotation occurs naturally 
over the course of the year and distributed the obser- 
vations as uniformly over the observing cycle as the 
scheduling allowed. Given the observing windows for 
Fomalhaut, this led to observations being acquired in 
August and December of 2010, and January and July 
of 20 1 1 , with a range of position angles as summarized 
in Table Q] 



Basic data reduction for all the observations was 
performed with the Spitzer Science Center (SSC) 
IRAC Pipeline (version SI 8. 18.0), which produced 
Basic Calibrated Data (BCD) frames and data qual- 
ity masks for each individual exposure. We also 
used the p ost-BCD IRACproc package (version 4.3; 



Schust er et al . 2006) for the sole purpose of remov- 



ing outliers (cosmic-rays) for each frame. Subsequent 
steps were performed with custom procedures in IDL. 
An extra bad pixel removal step was introduced in 
order to identify and remove residual bad pixels that 
occurred only in single frames. This was done by iden- 
tifying outliers from the median of each quadruplet of 
frames that were taken contiguously for a given dither 
position. The absolute center of the PSF in each frame 
was determined by cross-correlating the spider pattern 
with itself after a rotation by 180 degrees. All frames 
were then shifted to a common center and oversam- 
pled to a pixel scale of 300 mas/pixel, after which the 
ADI-based PSF subtraction could commence. This 
procedure is based on the Locally O ptimized Combi 



nation of I mages (LOCI) procedure dLafreniere et al 



2007 , 20091) . with adaptations to suit the new type of 



observational scheme and the science goals. The main 
aspects of the data set that distinguish it from most 
types of situations in which LOCI is applied are the 
extremely high PSF stability for space-based observa- 
tions, the very large number of frames (48 x 8 = 384), 
and the relatively small separation of the region of 
primary interest (compared to the PSF size, FWHM of 
1 .72"). Particularly the latter two factors make the data 
prone to substantial self-subtraction in a regular LOCI 
context. Since one of our main objectives in this study 
is to set a firm upper flux limit in case no detection 
would be made, we adapted the procedure in such a 
way as to avoid self-subtraction to a very high degree, 
while still maintaining a strong contrast performance. 
This is a conservative approach, and we note that it is 
entirely possible that the contrast could be even further 
enhanced with a more aggressive LOCI implementa- 
tion. This will be a subject of future studies. 

Our adapted LOCI implementation follows the fol- 
lowing procedure: First, each frame is sequentially 
subjected to an individual optimized reference PSF 
construction and subtraction. The basic optimization 
area is an annulus with an inner radius of 15 pixels 
and an outer radius of 60 pixels, centered on the star. 
This optimization area is split up into pieces of approx- 
imately 9-by-9 pixels (~1.5 FWHM). In practice, this 
is done in such a way that the basic annulus is split 
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up into four annuli, each 9 pixels in width, and each 
annulus is split up azimuthally in such a way that an 
integer number of segments of equal angular width are 
created, and such that the length of the inner edge of 
the segment is as close to 9 pixels as possible. We 
will refer to these segments as 'exclusion areas' hence- 
forth. For each area, an individual LOCI optimization 
is executed where the optimization area consists of the 
full 15-60 pixel annulus, but with the exclusion area 
removed. Based on the resulting LOCI coefficients, 
a subtraction is then made only in the exclusion area. 
The resulting small area is saved to a frame which is 
put together piece by piece from the subtractions cor- 
responding to the respective exclusion areas. Hence, in 
this procedure, the subtraction area is equal to the ex- 
clusion area, and is completely non-overlapping with 
the optimization area. The advantage of this approach 
is that it becomes impossible for the algorithm to sys- 
tematically fit for any companion in the data, and thus 
the self-subtraction will be approximately zero. The 
possible cost comes from the fact that we exclude the 
stellar PSF regions that are physically the closest to 
the companion, and which may therefore correlate best 
with the actual PSF noise at the exact position of the 
companion. As we will see in the testing described in 
Sect. [3] the procedure indeed works extremely well for 
avoiding any self-subtraction. 

For the procedure described above, the reference 
frames are chosen from the available library of frames 
based on how far separated the position angles are 
between the subtraction frame and a given reference 
frame. Since we wish to ensure a very low degree 
of self-subtraction of any real companion, we conser- 
vatively choose that the separation must be at least 
1.72"(1 FWHM) at the inner edge of the exclusion 
area. This corresponds to different position angles at 
different separations from the star, hence the subtrac- 
tion areas at larger separation typically have access to a 
larger number of reference frames. However, since the 
position angles are spread over a large range (see Table 
Q3, every examined position has access to a sufficient 
number of reference frames for a very high quality PSF 
subtraction. 

Once all frames have been subjected to the PSF sub- 
traction, they are de-rotated to a common sky orien- 
tation where North is up and East is to the left, and 
collapsed into a final frame using the median of the 
individual frames. Another step is then performed, 
in order to take boundary effects into account. If a 
real companion happens to be positioned exactly on 



the boundary between two or more exclusion zones, it 
is partially vulnerable to residual oversubtraction even 
in the above procedure. To overcome this effect, we 
perform three additional LOCI procedures in exactly 
the same way as described above, with only one dif- 
ference: In the first additional procedure, the exclu- 
sion zones are shifted one half step in the azimuthal 
direction, in the second one, they are shifted one half 
step in the radial direction, and in the final one, they 
are shifted half a step in each direction. In this way, 
four different reduced frames are available to check 
if some companion gets fainter due to boundary ef- 
fects in some of the frames, and additionally allows to 
check for spurious features that could occur in some 
reductions but not in others. In this case there is a box- 
like artefact in two of the images but otherwise they all 
show nicely consistent patterns, hence for the further 
analysis we use the mean of the frames from the four 
reductions. 

In order to evaluate the achieved contrast as a func- 
tion of se paration, we use the sa me procedure as de- 
scribed in Marengo et al. ( 2009b . evaluated from the 
standard deviation in consecutive 1 -pixel annuli. An 
important difference is that we always use a 5a cri- 
ter ion here for our meas urements, rather than 3a as 
in iMarengo et al. d2009h . This is more stringent in 



the presence of speckle noise, although note also that 
there is a fairly close equivalence between a 5a single- 
signature criterion and a 3a double-signature criterion 
( Janson et al. 2008), the latter of which is relevant for 
the Mareng oetalT(l2009l) data. 



3. Results and Analysis 

We show the final reduced image in Fig. Q] and 
the corresponding sensitivity limits in Fig. where 
we also plot the expected brightnesses for planets of a 
few different masses at ages of 200 and 400 Myr. A 
very substantial improvement in contrast performance 
is achieved with our LOCI implementation, with an 
order of magnitude improvement in flux detectabil- 
ity compared t o the conventional ADI reduction in 
Mar engo et al. (2009), and with a more stringent de- 
tection criterion. No signature is found at the position 
where Fomalhaut b would be expected. Hence, we 
estimate an upper limit (5a) at the relevant position 
based on the standard deviation in a 7-by-7 pixel box 
(~1.2 FWHM on the side). This gives an upper limit of 
16.7 mag, which corresponds to 38.8 /iJy, again more 
than an order of magnitude improvement over previous 
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Table 1: Log of Fomalhaut observations in Spitzer program 70009. 



AOR ID 


Chan. 


Exp. 


Frames 


MJD 


PA 


40250112 


4.5 


10.4 s 


48 


55416.265 


254.5° 


40249856 


4.5 /im 


10.4 s 


48 


55423.457 


257.6° 


40249600 


4.5 /Ltm 


10.4 s 


48 


55547.666 


52.6° 


40249344 


4.5 


10.4 s 


48 


55561.778 


58.4° 


40249088 


4.5 /im 


10.4 s 


48 


55571.925 


62.2° 


40248832 


4.5 /urn 


10.4 s 


48 


55579.832 


65.0° 


40248576 


4.5 /Ltm 


10.4 s 


48 


55762.480 


244.5° 


40250368 


4.5 ^m 


10.4 s 


48 


55765.605 


245.7° 




Fig. 1 . — Final reduced image for the real data (left) and for the data with an artificial companion introduced at the 
expected position of Fomalhaut b (right). Arrow 1 points out the expected position of the companion based on earlier 
detections in the visible-light images. There is no corresponding point source seen in the real data. The artificial 
companion in the right-hand image passes through the data reduction without any flux loss, verifying that the non- 
detection is real, such that a stringent upper flux limit can be set. Arrow 2 points toward the brightest possible point 
source in the field. Its position is consistent with a ring-nested orbit, but the significance is too low to make any 
assessment of whether or not it is a real object. North is up and East is to the left in the images. 



4 



15 


'2 Mjup, 200 Myr 




15.5 






D) 
CO 

£ 16 

>< 

> 


2 Mjup, 400 Myr 
1 Mjup, 200 Myr 




w 

S 16.5 

CO 

e 


1 Mjup, 400 Myr 




17 


' 0.5 Mjup, 200 Myr 




17.5 


' 0.5 Mjup, 400 Myr 





6 7 8 9 10 11 12 13 14 

Separation (arcsec) 



Fig. 2. — Sensitivity limit as function of separation from Fomalhaut. The solid line is the azimuthally averaged sensi- 
tivity profile. There is a bump around 12"due to concentrated noise at some ranges of position angles. The expected 
position of Fomalhaut b from the visible wavelength range detections is in a cleaner part of the image space, and the 
local sensitivity at this position is shown by the black asterisk. Also plotted as dashed lines are t he expected bright- 
nesses for planets with a few different combinations of mass and age, according to models from Spiegel & B urrows! 
d201ll) . 
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data. 

In order to test that the non-detection is real and 
not an effect of any unexpected oversubtraction in the 
LOCI procedure, we make a full reduction following 
the exact same procedure as described in the previ- 
ous section, except we also introduce a faint artificial 
companion (a Gaussian with 1.72"FWHM) in all the 
pre-LOCI frames, at the expected position of Fomal- 
haut b, with a flux of 57 /iJy (this corresponds to an 
effective temperature of ^250 K, or equivalently 0.5- 
1 M jup at 200 Myr and 1-2 M jup at 400 Myr). The 
artificial companion passes through the LOCI reduc- 
tion entirely unaffected, with the same measured flux 
in the final frame as the flux that was put in within the 
error bars, and is well visible at la confidence in the 
final frame (see Fig. Q]). Note that the introduction of 
a new feature in the data affects the LOCI reduction 
itself - since the artificial companion exists in most 
reference frames as an additional feature that has to 
be fit for by the algorithm, the fit quality will in gen- 
eral be slightly worse on average. Since the extent of 
the companion is small with respect to the optimiza- 
tion area and it is faint with respect to the stellar PSF, 
the effect is small but noticeable as a marginally higher 
general noise level in the final reduced frame. In sum- 
mary, the procedure described here validates that the 
non-detection is real, and thus the upper flux limit is 
relevant. 

We show our upper limit in Fig. [3j along with 
other upper limits from the literature at various wave- 
lengths, as well as the detection values in the visual 
wavelength range. We also compare these values with 
vario us theoretical models . One model spectrum is 
from lBurrows et all d2003l) (henceforth BSL03), cor- 
responding to a 2-3 Mj U p planet (interpolated be- 
tween 2 Mj Up and 5 Mj up to match the flux value 
at F814W3) a t ~20 Myr. As noted in K08 and 
Marengo et al. 2009), his model is similar to e.g. the 
Fortnev et al.1 d2008l) model except in H-band where 
the BSL03 flux is higher. It can be seen from the 
comparison to the data that in the context of this 
model, a thermal flux interpretation is entirely incon- 
sistent with both the non-detection in H-band and our 
non-detection in the IRAC 4.5 /im channel. The H- 
band flux is strongly model-dependent as it is sensi- 
tive to uncertainties in opacity and the treatment of 



clouds. However, by contrast, the 4.5 fim flux is 
very insensitive to these effects, and varies only very 
marginally across differe nt models. To show this , 
we use n ewer models fro m Spiegel & Burrows! d201ll) 
based on iBurrows et al.1 d201 lh (henceforth BHN11). 
which has improved opacities and cloud treatment, and 
also includes a set of entirely cloud-free models for 
comparison. 

By using these new models and including clouds, it 
is possible to suppress the H-band flux to a significant 
extent (though typically still not quite to a sufficient 
extent for a non-detection to be consistent with ther- 
mal flux in F814W). However, this is not the case at 
4.5 fim. The flux remains very stable for a constant ef- 
fective temperature (typically ~400 K for these model 
comparisons), regardless of how the clouds are treated 
including the cloud-free case, independently of opacity 
treatment and also of specific metallicity (the BHN11 
models provide both Solar and super-Solar metallic- 
ity cases). We conclude that the effective temperatures 
required to get any substantial contribution of thermal 
flux to the observed F814W data point are simply in- 
consistent with the non-detection at 4.5 /im. 

In order to comply with the upper limit at 4.5 /im, 
we need effective temperatures in the range of ^200 K 
or lower, corresponding to, e.g., 1 Afj up at 400 Myr 
(see Fig. [4j. As a side point, this latter age would 
correspond to a new but so far unpublished estimate 
(E. Mamajek, priv. comm.) which is a bit older 
than the mean estima te of 200 Myr used in K08 and 
Mar engo et al. (2009). Here we do not make any as- 



1 F606W and F814W are filters for the Advanced Camera for Surveys 
on the Hubble Space Telescope, centered on wavelengths of ~600 
and ~8()0 nm, respectively 



sessment of the relative credibility of these two es- 
timates, but simply remark that it has no real rele- 
vance for the spectral comparisons we are perform- 
ing here. The main factor (beyond cloud treatment 
and opacity) that affects the spectral energy distribu- 
tion is the effective temperature. The mass to which 
this temperature corresponds depends on the age and 
vice versa, but this has little impact on the spectrum, 
especially for the small discrepancy in age that we are 
concerned with here. Hence, as long as we do not ac- 
tually try to determine the mass, we do not need to 
assess which of 200 Myr or 400 Myr is the better esti- 
mate. While on this note, it is worth pointing out that 
hot/warm/cold-start models are of no significant rele- 
vance to this discussion, since convergence will have 
occurred at these ages and masses, regardless of initial 
entropy dSpiegel & Burrowsll201 1 ). 

For the colder models required to match the 4.5 /jm 
upper flux limit, virtually all flux at shorter wave- 
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Fig. 3. — Model comparisons to the observational data for a ^400 K atmosphere (~2-3 Mj up at ~200 Myr for 
BSL03, 4 Mj U p at 200 Myr for BHN11). The solid line in each panel is the model spectrum, and the black lines 
are the corresponding fluxes in the respective photometric bands. Crosses mark detections and triangles mark upper 
limits. HST data points from K08 are in green and Keck/Gemini data points from K08 are in brown. Note that the 
error bars at F6 06W and F814W are smaller than the sizes of the symbols. Red symbols are Spitzer upper limits from 
Marengo et al. (2009). Our new upper limit from Spitzer is the magenta triangle. Upper left: A BSL03 model. Upper 
right: A BHN11 model with patchy clouds and Solar composition. Lower left: A BHN11 model with patchy clouds 
and increased metallicity. Lower right: A BHN1 1 model with Solar composition and a cloud-free atmosphere. Models 
in this effective temperature range are required to produce an adequate amount of flux in F814W, but they are typically 
inconsistent with the upper flux limits in H-band and/or L'-band, and always fully inconsistent with the upper limit at 
4.5 fim 
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Fig. 4. — Model comparisons to the observational data for a ^200 K atmosphere (1 Mj up , 400 Myr). The symbols 
have the same meaning as in Fig. [3] Upper left: A BSL03 model. Upper right: A BHN1 1 model with patchy clouds 
and Solar composition. Lower left: A BHN11 model with patchy clouds and increased metallicity. Lower right: A 
BHN1 1 model with Solar composition and a cloud-free atmosphere. These models are marginally consistent with the 
upper flux limit at 4.5/im, but predict one to several orders of magnitude too little flux to be consistent with thermal 
radiation in F814W. 
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lengths is lost, and there is no way to match the F8 14W 
point with any thermal flux. The closest case is the 
extreme case of a completely cloud-free atmosphere 
and solar metallicity, but also in this case the flux is 
more than an order of magnitude too small (see Fig. 
@). Increasing metallicity has the effect of decreasing 
flux at 800 nm compared to 4.5 /im, hence decreasing 
metallicity could have the opposite effect. However, 
Fomalhaut A has a super-solar metallicity with a mean 
measured valu e of 0.3 dex in a compilation of litera- 
ture values in Soubiran d2010h . Fomalhaut b should 
have an equal or larger metallicity - thus, adjusting the 
metallicity is not a feasible route toward reaching con- 
sistency with the observational data. 

Based on the above results and considering the fur- 
ther aspects of the collected body of observations of 
Fomalhaut as will be discussed in detail in Sect. [4] 
it is highly unlikely that the observed flux at visible 
wavelengths has any direct connection to the suspected 
giant planet that might shepherd the debris disk of Fo 



malhaut and force it into an eccentric state ( Kalasetal 



2005). In this context, and considering our very strong 
detection limits, it is interesting to assess whether 
this shepherding planet (what might be referred to as 
the 'real' Fomalhaut b) can be seen in our images. 
Since the shepherding planet can be as low in mass as 
0.5 Mj up dChiang et al] |2009). it is fully plausible that 
it could remain undetected in our images if the age is 
as old as 400 Myr, and since its orbit covers a range 
of projected separations, it can also hide in some parts 
of the orbit even if the mass is slightly higher. How- 
ever, we do cover a very large fraction of its possible 
parameter space, and it is noteworthy that the brightest 
possible point-source in the field with a significance 
of 4.3ct is in fact located at a position that would be 
consistent with a ring-nested orbit (arrow 2 in Fig. Q]). 
However, the fact that the possible point source is not 
at a 5 a confidence level obviously means that more 
data or an even further improved PSF reduction would 
be necessary in order to test its validity. 

4. Discussion 

In this section, we discuss in detail how our non- 
detection at 4.5 /im affects the interpretation of Fo- 
malhaut b, and what can be deduced about the system 
from the full body of existing data. The data points that 
exist for the detected point-source in K08 are two de- 
tections in the F606W filter from 2004 and 2006, and 
one detection in the F814W filter from 2006. These 



data points are shown in Figs. [3] and [4] Note that the 
object is variable between 2004 and 2006. This is not 
a small effect; in fact, it corresponds to a change in 
brightness (dimming) of a factor 2, at a confidence of 
8ct. This is the same level of confidence as, e.g., the 
second-epoch F606W detection of the point-source al- 
together. Hence, to the extent that we can trust the 
data at all, we must consider this variability as a real 
effect, and it needs to be accounted for in a comprehen- 
sive interpretation of the object. In addition to these 
data, a third epoch HST observation has also been ac- 
quired but has not yet been published at the time of 
writing (P. Kalas, priv. comm.). There are also up- 
pe r flux limits from n on-detections in a range of bands 



m lKalas et alj d2008l) : F435W, H, CH4S, CH4L, and 
L', and upper limits at 3.6, 4 .5, 5.8 and 8.0 /im from 
Spitzer ( Marengo et alj[2009b . the second of which we 
improve on in this article. 

There are two main lines of interpretation of the 
point source in K08, only one of which actually in- 
cludes any flux from a planet. In this scenario, the flux 
at F814W originates from the thermal emission of a 
planet, which has to be close to a mass of ^3 Mj up in 
order to fit the data point assuming an age o f 200 Myr 
(e.g. lBurrows et al.ll2.Q03h iFortnev et al.ll2.008b . This is 
however poorly consistent with the rest of the avail- 
able data. Most obviously, it is a factor 20^-0 brighter 
in F606W than expected in such a scenario. In order 
to explain the F606W data in terms of brightness and 
variability, K08 infer a hypothesis of Ha accretion on 
the planet. Given that this would require gas accre- 
tion at the same rate of a few Myr old TTauri stars like 
GQ Lup (accretion rate calculated in the supplemental 
material of K08), whereas Fomalhaut is 200^-00 Myr 
and has no other known signs of gas anywhere in the 
system, we consider this hypothesis highly unlikely. 

Aside from this, there is the issue that a ^3 Mj up 
companion should also be detectable in the near- 
infrared. This was an issue already in K08, where 
the H-band flux predicted by theoretical models was 
significantly higher than the upper limit from obser- 
vations, and became an even larger is sue given the 
Spitzer data published bv lMarengo et al.l (120091) . where 
the additional limits at larger wavelengths provided 
very little opportunity for any flux to remain unde- 
tected from such a companion. As was shown in the 
previous section, our new Spitzer data now provides 
even much tighter constraints on the thermal emission 
hypothesis, and provides the opportunity to conclu- 
sively address this issue. We find that any thermally 
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emitting companion responsible for the F814W flux 
(effective temperatures of ^400 K, e.g. 2-3 Afj up 
at 200 Myr) would have emitted more than an order 
of magnitude more flux at 4.5 /im than our 5<7 upper 
limit, regardless of the choice of theoretical models. 
Conversely, any companion that thermally emits ra- 
diation at 4.5 /zm at levels comparable to our upper 
limit (effective temperatures of ^200 K, e.g. 1 Mj up 
at 400 Myr) would emit at least an order of magni- 
tude too little flux at F814W to explain the observa- 
tions, and for most realistic model parameters (e.g. 
any inclusion of clouds) the flux would be even much 
smaller. Hence, we can firmly exclude the hypothesis 
that any of the observed flux in K08 actually originates 
from a giant planet. 

Given that the SED of the K08 point source can be 
interpreted as having a reasonable match to the stel- 
lar SED, it seems more likely that what is seen in the 
K08 images is some form of reflected or scattered ra- 
diation from the star. Given the large effective area 
that is required for this, the only plausible origin for 
reflected/scattered radiation is dust. Hence, it is likely 
that we are seeing a concentration of dust, which may 
or may not be associated with a planet. We will dis- 
cuss some dust-related interpretations in the following: 
firstly, we will consider the hypothesis of an optically 
thick disk around a giant planet. This is the second pre- 
ferred scenario in K08, and requires a disk radius of at 
least ^20 i?j U p (for the case of high albedo; larger for 
the case of low albedo). It may not be unreasonable 
that circumplanetary rings start out with such proper- 
ties. However, it might also be argued that the age 
of the Fomalhaut system of 200^1-00 Myr should have 
provided adequate time for moons to form and exca- 
vate the disk, leaving only rings of a much smaller ef- 
fective size. Regardless of whether the optically thick 
disk scenario is fundamentally realistic or not, there 
are several reasons for why this scenario is inconsistent 
with the existing data. One very important constraint 
is that such a disk cannot account for the factor 2 vari- 
ability observed in F606W. Furthermore, the spin of 
the st ar has been recently m easured with interferome- 
try dLe Bouquin et al . 2009). If the spin of the star is 
aligned with the plane of the disk, this means that the 
fainter Western side is closer to us than the brighter 
Eastern side. Although this is opposite to what would 
be expect ed for purely for ward-scattering dust, it is 
shown in iMin et al.1 (2010) to be consistent with the 
scattering behaviour of large dust grains. If it is indeed 
the case that the Western side of the disk is closer to us, 



and the K08 point source orbits within the disk, then 
it follows that the object is located between the par- 
ent star and Earth, in the radial direction. It would be 
very difficult for an optically thick disk to reflect large 
amounts of light to Earth under such circumstances. 

In addition to these points, another important argu- 
ment against the involvement of a giant planet (which 
also applies to the thermal emission case discussed 
above) is the orbit of the object. The third epoch astro- 
metric measurement implies a ring-crossing orbit for 
the point source (P. Kalas, priv. comm.). Although this 
is not yet published, it would not be surprising given 
the previously published data, because, in fact, already 
the first two epochs are inconsistent with an orbit that 
traces the edge of the ring, at a ^2er level (the direction 
is largely consistent with such an orbit, bu t the sp eed is 
not). This is pointed out by lChiang et al. (2009), who 
do not put a large emphasis on this fact, as they argue 
that the error bars are probably underestimated. How- 
ever, we note that this is diametrically opposite to the 
interpretation of the astrometric errors in K08, where 
it is concluded that they are an upper limit to the real 
error (K08, supplemental material). A ring-crossing 
orbit would be inconsistent with an association of the 
point source to a giant planet, as it would strongly af- 
fect the geometry of the disk, hence any planet asso- 
ciated with the point source would need to be low in 
mass. 

In the context of a low-mass planet, we note that 
one way to suppress the 4.5 /zm flux with respect to 
shorter wavelengths is to consider hotter temperatures 
and smaller surfaces. Effective temperatures higher 
than the ^200^400 K that we have considered thus far 
are not reasonable for isolated objects at the system 
age, particularly for smaller (and thus lower-mass) ob- 
jects than Jupiter-class planets. However, following 
intense bombardment of planetesimals, rocky proto- 
planets may acquire molten surfaces and reach effec- 
tive temp eratures of ~ 1000-300 K over brief periods 
of time (Miller-Ricci et al . 2009). Hence, the possibil- 
ity that the observed light-source could be a hot colli- 
sional afterglow of a < 10 MEarth object should not be 
dismissed out of hand. Still, it would probably be dif- 
ficult to reproduce all the observed data points in such 
an interpretation, particularly the detection in F606W, 
and the simulta neous non-detections in H -band and 
at 4.5 /im. The iMiller-Ricci et al. (2009) models do 
not cover the full relevant wavelength range, but if 
we consider, for example, a planet of 10 Mgarth an d 
1.8 i?Earth and work from pure blackbody consider- 
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ations, we can establish that the brightness temper- 
ature required to reproduce even the lowest of the 
F606W data points (0.30 /Jy) is more than 1500 K, 
whereas the upper limits in H-band (0.71 /iJy) and at 
4.5 fim (38.8 /iJy) both require brightness tempera- 
tures of ^700 K or less. We consider that future mod- 
elling efforts would be worthwhile to examine whether 
these conditions and the rest of the flux limits can all be 
simultaneously fulfilled, but for the purpose of our dis- 
cussion here, we simply treat it as an option that can- 
not be categorically excluded. On balance, it should 
be noted that an observation of this type of scenario 
is probably rather unlikely, for several reasons, includ- 
ing the fact that it is expected to last over timescales of 
^10 4 yr for the case of a thin atmosphere, very short 
compared to the age of Fomalhaut. The timescale can 
be extended if the atmosphere is thickened and clouds 
are included, but the observable brightness tempera- 
ture decreases accordingly. As a side note, if clouds 
were involved, they could possibly account for the 
variability in F606W in this scenario. 

Given the SED of the point source and its vari- 
ability, along with the considerations above, the per- 
haps most plausible way to consistently explain the ob- 
served properties of the K08 point source is through a 
cloud of dust, which is either transient or has a tran- 
sient component. There are two possible scenarios 
associated with such an interpretation that have been 
suggested in the literature. In one scenario, the ob- 
served point source is a residual (gradually dispersing) 
dust cloud from a recent planetesimal collision. We 
certainly know that such collisions should occur fre- 
quently in the Fomalhaut system, given that they are 
the very origin of dust in debris disks. This scenario is 
mentioned by K08, who argue against it based on the 
fact that such collisions should be much more com- 
mon within the actual ring feature than just outside of 
it where the point source is observed, hence the relative 
probability to observe it where it is observed should be 
low. This is certainly true, but we note that there is 
a clear selection effect involved - due to the high vi- 
sual brightness of the ring and the speckle-like nature 
of the noise, any number of equivalent events that hy- 
pothetically do happen within the ring feature would 
be likely to pass unnoticed. One might also hypothet- 
ically imagine that the cloud is in the present position 
for some specific dynamical reason, for instance if the 
material is trapped in resonance with a giant planet sit- 
uated elsewhere. 

The second scenario is essentially the same as the 



first, but involves a central rocky /icy object with a mass 
less than ~10 Msarth, to which a swarm of plan- 
etesim als is gravitationally bound jKennedy & Wyatt 
|201 lb . Collisions between these planetesimals pro- 
duce the observed dust. We consider both of these sce- 
narios to be reasonable within the constraints set by the 
data, and simply conclude that the K08 point source is 
well consistent with a transient or semi-transient dust 
cloud, which may or may not be gravitationally bound 
to a central object of planetary mass. With regards to 
the fact that the point source has been frequently re- 
ferred to as a directly imaged planet in the literature, 
we note that this is incompatible with the observational 
evidence, for two independent reasons: ( 1 ) Although it 
cannot be formally excluded, there is insufficient evi- 
dence to support that there is any compact object of 
planetary mass associated with the point source alto- 
gether. (2) Even if such an object is present, in sev- 
eral of our considered scenarios we do not observe any 
photons from this object itself, hence it cannot be es- 
tablished that it has been directly imaged. 

5. Conclusions 

In this paper, we have presented observations per- 
formed with SpitzerflRAC in the 4.5 /um band for 
the purpose of trying to detect thermal emission from 
Fomalhaut b. A new LOCI-based PSF subtraction 
scheme was implemented to achieve high contrast with 
minimal companion flux loss, which enabled an or- 
der of magnitude improvement in contrast-limited sen- 
sitivity with respect to previous efforts. The non- 
detection of any flux at the expected position can there- 
fore be used to provide strong constraints on the un- 
derlying physics of the point-source seen at visible 
wavelengths. In particular, we find that there is al- 
most certainly no direct flux from a planet contribut- 
ing to the visible-light signature. This, in combination 
with the existing body of data for the Fomalhaut sys- 
tem, strongly implies that the dynamically inferred gi- 
ant planet companion and the visible-light point source 
are physically unrelated. This in turn implies that the 
'real' Fomalhaut b still hides in the system. Although 
we do find a tentative point source in our images that 
could in principle correspond to this object, its signifi- 
cance is too low to distinguish whether it is real or not 
at this point. 

Concerning the visible-light point source, its under- 
lying physics is unclear, but the only hypothesis that 
can be shown to reasonably fit all existing data is an 
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optically thin dust cloud, which is transient or has a 
transient component. If this interpretation is valid, the 
cloud may or may not be physically bound to a central 
object in the super-Earth mass regime. 

The Fomalhaut system is a rich topic of conver- 
sation, and we thank Adam Burrows, Carsten Do- 
minik, James Graham, Ray Jayawardhana, Paul Kalas, 
Michiel Min, and many others for interesting discus- 
sions. This work is based on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. M.J. is funded by 
the Hubble fellowship. J.C.C., J.R.B., and PW. were 
supported by grant AST- 1009203 from the National 
Science Foundation. D.S.S. gratefully acknowledges 
support from NSF grant AST-0807444 and the Keck 
Fellowship. 
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